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Abstract 

We have performed high precision measurements of the zero-energy neutron scattering amplitudes of gas phase 
molecular hydrogen, deuterium, and 3 He using neutron interferometry. We find b np — (—3.7384 ± 0.0020) fm[l], 
b nd = (6.6649 ± 0.0040) fm[2,l], and 6 n3 He = (5-8572 ± 0.0072) fm[3]. When combined with the previous world 
data, properly corrected for small multiple scattering, radiative corrections, and local field effects from the theory of 
neutron optics and combined by the prescriptions of the Particle Data Group, the zero-energy scattering amplitudes 
are: b np = (-3.7389 ± O.OOfO) fm, b nd = (6.6683 ± 0.0030) fm, and b n3Ue = (5.853 ± .007) fm. The precision of 
these measurements is now high enough to severely constrain NN few-body models. The n-d and n- 3 He coherent 
neutron scattering amplitudes are both now in disagreement with the best current theories. The new values can be 
used as input for precision calculations of few body processes. This precision data is sensitive to small effects such 
as nuclear three-body forces, charge-symmetry breaking in the strong interaction, and residual electromagnetic 
effects not yet fully included in current models. 
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1. Introduction 

The last decade has seen a revolution in the accuracy 
with which low energy phenomena in nuclear few body 
systems can be calculated. Insight into certain features 
of few-nucleon systems has come both from greatly- 
improved calculations using potential models based 
on the measured nucleon-nucleon (NN) interaction [4] 
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and also from the development of effective field the- 
ory (EFT) approaches based on the chiral symmetry of 
QCD [5,6,7]. Such theories have been used to develop a 
physical understanding rooted ultimately in QCD for 
the relative sizes of many quantities in nuclear physics, 
such as nuclear N-body forces [8] and in particular the 
nuclear 3-body force (3N), which is now investigated in- 
tensively. Although it is well understood that 3N forces 
must exist with a weaker strength and shorter range 
than the NN force, little else is known. 

EFT has been used to solve the two and three nu- 
cleon problems with short-range interactions[7,9]. For 
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the two-body system, EFT is equivalent to effective 
range theory and reproduces its well-known results for 
NN forces [10,11,12]. The chiral EFT expansion does 
not require the introduction of an operator correspond- 
ing to a 3N force until next-to-next-to leading order in 
the expansion, and at this order it requires only two 
low energy constants [13,14], which are taken to be the 
triton binding energy and the zero energy doublet n-d 
scattering amplitude. There have also been significant 
advances in other approaches to the computation of the 
properties of few-body nuclei with modern potentials 
such as the AV18 potential [15,16], which includes elec- 
tromagnetic effects and terms to account for charge- 
independence breaking and charge symmetry breaking 
of the strong interaction. These calculations accurately 
reproduce the well-measured energy levels of few body 
bound states only with the phonomological inclusion 
of a nuclear 3-body force, so it is clear that more infor- 
mation on this force is needed for further progress. 

Precision measurements of low-energy strong inter- 
action properties, such as the zero-energy scattering 
amplitudes and electromagnetic properties of small A 
nuclei, are therefore becoming more important for low 
energy, strong interaction physics both as precise data 
that can be used to fix parameters in the EFT expan- 
sion and also as new targets for theoretical prediction. 
In this article, we briefly summarize a series of preci- 
sion measurements of the zero-energy coherent scatter- 
ing amplitudes in the two, three, and four-body sys- 
tems using neutron interferometric techniques[l,2,3]. 
We summarize these zero-energy coherent scattering 
amplitude results and discuss possibilities for future 
measurements. 



2. Neutron Optics Theory 

The zero-energy coherent scattering amplitude is 
the linear combination of scattering amplitudes that 
gives rise to the optical potential of a neutron in a 
medium[17]. The zero-energy bound scattering ampli- 
tude, b, is related to the free scattering amplitude a by 



and for 3 He, it is the linear combination of the singlet 
and triplet scattering amplitudes, 



m + M 



(1) 



Here, m is the mass of the neutron and M is the mass of 
the atom. For hydrogen, a is the linear combination of 
the singlet and the triplet scattering amplitudes given 
by, 



a„ p = (1/4) 1 a np + (3/4) 3 a n 



(2) 



for deuterium it is the linear combination of the doublet 
and quartet scattering amplitudes, 



«n3He = (V4) ^sHe + ( 3 /4) 3 a n 3 He - 



(4) 



The bound zero-energy scattering amplitude is one par- 
ticular linear combination of the triplet and singlet (or 
doublet and quartet) scattering amplitudes. Knowl- 
edge of some other combination allows one to indepen- 
dently extract the individual bound (or free) scattering 
amplitudes for each state. 

The phase shift measured in neutron interferometry 
is proportional to the real part of the S'-wave coherent 
scattering amplitude in the medium, 



. = fc(l-n)£>= —XNDb. 



(5) 



and = (1/3) a nd + (2/3) a n 



(3) 



where k is the incident wave vector, n is the index of 
refraction, N is the number density, D is the thickness 
of the sample, and A is the neutron wavelength. Thus 
to experimentally measure b, the neutron optical phase 
shift <f>, the atom density, the sample thickness, and 
the neutron wavelength must each be measured to high 
precision. 



3. Experimental Procedure 

Scattering amplitude measurements were performed 
at the National Institute of Standards and Technology 
(NIST) Center for Neutron Research (NCNR) Interfer- 
ometer and Optics Facility [18]. A cold monochromatic 
neutron beam (E = 11.1 meV) enters the perfect sili- 
con crystal neutron interferometer and is coherently di- 
vided via Bragg diffraction into two beams that travel 
along paths I and II as shown schematically in Fig. 1. 
These beams are again diffracted and then coherently 
recombined to form the interference pattern. A detailed 
description of the facility, experimental arrangement, 
and procedures for the determination of zero-energy 
neutron scattering amplitudes can be found in Ref. [1]. 

A secondary sampling method is used to measure 
the phase shift (j> due to the gas sample. This is accom- 
plished by positioning a rotatable quartz phase shifter 
across the two beams as shown in Fig. 1. The intensities 
of the beams that arrive at the two 3 He detectors are 
a function of the phase shifter angle S and are given by 

Io{S) =A + B cos(Cf(S) + 4> gaa + <t> cM ) (6) 

I H {S) = A H + B COs{Cf(8) + 4>gas + 4>cell + n). 

The values of Ao, Ah, B, and C are extracted from fits 
to the data. The function f(S) depends on the Bragg 
angle 6b and is a measure of the neutron optical path 
length difference between the beams induced by the 
phase shifter and is given by 
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Fig. 1. A schematic view of the experimental setup as the 
neutron beam passes through the perfect crystal silicon inter- 
ferometer. Parameters associated with the neutron optics are 
discussed in the text. 

sin(e B )sin(J-J ) . . 

/l j cos 2 (0 B )-sin 2 (<5-5 o )' 1 j 

The hydrogen, deuterium, or 3 He gases are housed in 
a cell specifically designed to minimize the phase shift 
4>ceii due to the aluminum walls of the cell (see Fig. 1). 
The phase shifts arising from the presence of the gas 
and cell (<fi g as and <j> C eii) were determined by collect- 
ing « 10 interferogram pairs with the cell positioned 
both within the interferometer and removed from the 
beam paths. The phase difference between cell-in/cell- 
out sets of interferograms is extracted for each pair, 
with a typical set shown in Fig. 2. The phase shift from 
the cell, (ficeii, was determined using an evacuated cell. 

The atom density was determined using the mea- 
sured purity of the gas and the ideal gas law with 
virial coefficient corrections up to the third pressure co- 
efficient. The absolute temperature was continuously 
monitored using two calibrated 100 fl platinum ther- 
mometers that have an absolute accuracy of 0.023 % 
at 300 K. The pressure was continuously monitored 
using a calibrated silicon pressure transducer capa- 
ble of measuring the absolute pressure to better than 
0.01 %. The wavelength of neutrons traversing the in- 
terferometer was measured using a pyrolytic graphite 
(PG 002) crystal. This analyzer crystal, calibrated sep- 
arately against a Si crystal with a precisely-known lat- 
tice constant, was placed in the H-beam of the inter- 
ferometer and rotated so that both the symmetric and 
anti-symmetric Bragg reflections were determined. In 
a separate test, the stability of the wavelength over the 
measurement time was shown to be 0.001 %. More de- 
tails on the measurement techniques and systematic 
uncertainties involved in the determination of the neu- 



Fig. 2. A typical pair of interferograms with 3 He present in 
the cell. The oscillations arise from the change in path lengths 
created as the phase shifter is rotated (see Fig. 1). Data are 
shown for both the cell in and out of the interferometer. 

tron wavelength, atom density, temperature, and cell 
thickness can be found in refs. [1] and [3]. 

The value of the bound zero-energy scattering am- 
plitude was calculated for each data set on a run- 
by-run basis. These values were combined using a 
weighted average to obtain b for each gas species. Our 
reported results are b np = (-3.7384 ± 0.0020) fm[l], 
b nd = (6.6649 ± 0.0040) fm[2,l], and & n3He = 
(5.8572 ± 0.0072) fm[3]. 



4. Results and Discussion 

Since these measurements of the n-p, n-d, and n- 3 He 
zero-energy scattering amplitudes were performed in 
an identical manner using the same apparatus (cell, 
neutron wavelength analyzer, and pressure and tem- 
perature monitors) , one can take the ratio of the mea- 
surement values to obtain results which are even less 
sensitive to any potential remaining systematic errors. 
Using b np as the reference, we obtain the ratios 

b n3Re /b np = (-1.5668 ± 0.0021) 
bnd/b np = (-1.7828 ±0.0014). 

Although these ratios possesses slightly larger statis- 
tical uncertainties, it will be independent of any un- 
known systematic uncertainty to first order and can 
provide a more robust target for comparison to the- 
ories which attempt to calculate all of the scattering 
lengths within a common theoretical framework. 

Recently, theoretical predictions of the n-d zero- 
energy scattering amplitude were performed using the 
high precision NN forces CD Bonn 2000, AV18, Nijm 
I, II and 93 in combination with 3N force models[19]. 
The results of these calculations are summarized in 
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Fig. 3. Theoretically determined values for the n-d scattering 
amplitude b nd and the triton binding energy using different 
NN and 3N forces[19]. The experimental value consists of our 
measurement of b nd and the very precisely known value of the 
triton binding energy, (8.481855 ± 0.000013) MeV[20]. 
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Fig. 4. Comparisons of the n- 3 He coherent and incoherent scat- 
tering amplitudes as determined from our measurement and the 
measuremetn of Zimmcr et al. [22] to theoretical calculations 
of these same parameters using either an R-matrix formalism 
or the AV18 NN potential with two types of 3N forces[21]. 



Fig. 3 alongside the experimental data. For NN forces 
alone with and without electromagnetic interactions, 
they recovered the approximate correlation between 
the triton binding energy and 2 a n d known as the 
Phillips line. The approximate correlation between 
these observables is now understood to be a generic 
feature for systems such as the deuteron and other 
low A nuclei whose size is significantly larger than the 
f fm scale of the NN interaction range. However it is 
understood that this correlation is only an approxi- 
mation and should not be obeyed exactly: for example 
3N forces will introduce corrections. Although the ad- 
dition of 3N forces of a wide variety of types does shift 
the values closer to the observed 2 a nd , none of these 
calculations agrees with the new high precision mea- 
surements. The authors note that the 2 a nc i zero-energy 
scattering amplitude must be considered as an inde- 
pendent low energy observable for future calculations 
in few body systems. 

A second group has recently published new calcu- 
lations of the spin-dependent n— 3 He scattering am- 
plitudes using the Resonating Group Method and a 
variety of modern NN and 3N potentials[21]. Com- 
parisons of the experimental free nuclear singlet and 
triplet scattering amplitudes as determined from our 
measurement and the n- 3 He incoherent scattering am- 
plitude from Zimmer et al. [22] with theoretical calcu- 
lations of these parameters were performed using an 
R-matrix formalism and the AV18 NN potential with 
two types of 3N forces and shown in Fig. 4. It is clear 
that the theoretical predictions from even these im- 
pressive, state-of-the-art calculations using resonating 
group techniques still lie outside the range of experi- 
mental uncertainties for the 3 He binding energy and 
the singlet and triplet n- He scattering amplitudes. 



The precision of these measurements is now high 
enough to severely constrain these few-body models. 
Both the n-d and n- 3 He coherent neutron scattering 
amplitudes are in disagreement with the best current 
theories. Only models which correctly take into account 
nuclear three-body forces, charge-symmetry breaking, 
and residual electromagnetic effects have the possibil- 
ity to successfully confront the data. 

We can combine these measurements with the previ- 
ous world's data to obtain new values for the coherent 
scattering amplitudes . We note that some of the past 
measurements of high precision using a gravity reflec- 
tometer performed in Garching a few decades ago [23] 
must be corrected for local field effects[24], which are 
significant for the coherent neutron scatting amplitude 
in hydrogen b np as measured in the gravity refrac- 
tometer but are negligible in neutron interferometry. 
This correction shifts the value from b np = (—3.7409 ± 
0.0011) fm[23] to b np = (-3.7390 ±0.0011) fm[24], an 
approximately two sigma effect. In addition there is 
another two sigma correction to the neutron interfer- 
ometry measurement in hydrogen gas due to multiple 
scattering and correlation corrections to the neutron 
index of refraction which are negligible for the gravity 
refractometer measurements. These corrections have 
already been applied to obtain the value of b np = 
(-3.7384±0.0020) fm[l] quoted above. After these cor- 
rections, both independent measurements are in excel- 
lent agreement. We can therefore present an improved 
value of the coherent neutron scattering length of hy- 
drogen of b np = (-3.7389 ± 0.0010) fm. These correc- 
tions turn out to be negligible for deuterium. 
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5. Future Measurements 

As mentioned above, the measured value of b np was 
corrected for multiple scattering and correlation cor- 
rections to the neutron index of refraction. A more pre- 
cise interferometric measurement in comparison with 
the gravity refractometer value would be able to isolate 
this term experimentally. Such a measurement would 
constitute the first experimental observation of cor- 
rections to the neutron index of refraction due to vir- 
tual excitations and multiple internal scattering within 
a molecule. Such effects were predicted by Nowak 20 
years ago[25]. The calculation of Nowak, which was 
done in the long wavelength limit kRo — > (Ro = bond 
length of H 2 = 0.74611 nm, for D 2 Ro = 0.74164 nm 
[26]), must be extended to the conditions of the ex- 
periment (which correspond to fci?o=T-73) to make a 
precise prediction. It would be possible to measure b np 
using H2 gas at least a factor of two more precisely, 
which may be accurate enough to isolate this neutron 
optics effect experimentally for the first time. 

Our measurement of b nd will hopefully soon be com- 
plemented by a measurement in progress at PSI of the 
incoherent n-d scattering amplitude using pseudomag- 
netic precession in a polarized deuterium target [27], 
which is also an interferometric technique that oper- 
ates in neutron spin space as opposed to real space. 
This measurement will allow one to separate the two 
spin channels and determine the interesting doublet 
amplitude, 2 a„d, with high precision. We have argued 
that the quartet amplitude should be independent of 
3N forces and the theoretical value should be reason- 
ably robust [28] and were thus able to extract a value 
of 2 a ni = (0.645 ±0.003(expt) ± 0.007(thcory)) fm[3], 
however an experimental determination is badly 
needed. We expect that these new measurements will 
improve 2 a nd by an order of magnitude to « 10~ 3 . 
The doublet amplitude is very important from a the- 
oretical point of view: in the EFT approach it fixes 
one of the two low energy constants in the expansion 
to next-to-next-to leading order (NNLO) and its de- 
termination will allow more precise predictions to be 
made for other few body systems. 

It is also interesting to consider how the experimen- 
tal determinations of the n- 3 He scattering amplitudes 
can be further improved. Our measurement of the zero- 
energy scattering amplitude is dominated by statis- 
tical uncertainties in the measurement of the phase 
shifts and consequently there is still room for improve- 
ment. In the case of the incoherent scattering ampli- 
tude (h) determination from pseudomagnetic preces- 
sion performed by Zimmer et al. [22], the accuracy is 
unfortunately limited by the poor experimental knowl- 
edge of the relative contributions of singlet and triplet 
channels to the n- 3 He absorption cross section. A bet- 



ter measurement of this ratio, currently known to ~ 
1 % [29,30], could be immediately combined with the 
Zimmer et al. measurement to improve the accuracy 
of bi by as much as a factor of three. In addition, a 
new measurement currently being planned at NIST to 
directly measure the spin-dependent n- 3 He scattering 
amplitudes is independent of this ratio [31]. Dramati- 
cally reduced uncertainties for *a 3jj e and 3 a n3 jr e in 
n- 3 He are therefore possible. 

Measurements in other light nuclei are also quite 
feasible. Two targets that we are presently exploring 
are 4 He and tritium. From a theoretical point of view 
at present, a 4 He measurement is not that interesting 
because the framework for solving five-body problems 
is not yet in place. The tritium measurement on the 
other hand is quite interesting theoretically because the 
small inelastic effects make it much easier to calculate 
than n- He, but it is experimentally difficult because 
of the radioactive nature of the target. We are investi- 
gating the possibility of performing a measurement of 
b nt using the same techniques used in the current mea- 
surements. We expect that the limiting factor in these 
measurements will be how well we can determine the 
isotopic composition of the radioactive tritium gas. 

Perhaps the single most interesting low energy 
NN scattering amplitude to measure is the neutron- 
neutron scattering amplitude b nn . A value for b nn is 
the last remaining obstacle to theoretically calculat- 
ing the n-A scattering amplitudes (A > 3) from first 
principles. At present, no direct measurements of b„ n 
exist. An experiment to determine b nn by viewing a 
high-density neutron gas near the core of a reactor and 
measuring a quadratic dependence of the neutron flu- 
ence on source power is currently being designed [32]. 
A second experiment to let the neutrons in an ex- 
tracted beam scatter from each other has also been 
considered [33]. An EFT analysis to extract b„ n from 
the iv~d — > nnj reaction has also recently been per- 
formed [34]. Any direct measurement would need to 
be performed to an accuracy of a few percent to be 
physically interesting. 

One additional measurement using neutron interfer- 
ometry is presently underway at NIST, a determina- 
tion of the neutron's mean square charge radius by 
measuring the neutron-electron scattering amplitude, 
&ne[35]. This experiment is expected to lead to a five- 
fold improvement in the precision of b ne . This precision 
may be sensitive to radiative corrections calculable in 
a recent EFT approach [36]. 

Theoretical advances over the last decade have made 
few-nucleon systems into a quantitative testing ground 
for low energy QCD. Low energy neutrons can be used 
to perform high-precision measurements of scattering 
amplitudes. We look forward to new high-precision 
measurements in this field in the next few years. 
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